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We consider a homogeneous non-ideal Bose gas at nonzero temperature in equilibrium below 
the critical temperature Tc in the framework of finite temperature field theory. An algorithm is 
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I. INTRODUCTION 



Experimental advances on Bose-Einstein condensation (BEC) in trapped gases of alkali atoms has caused a renewed 
interest in theoretical investigations of the phenomenon^ In particular, the weakly interacting Bose gas has received 
much attention, since it is interesting on its own right and serves as a simple model systerrP. 

While the basic properties of BECs at zero temperature are fairly well understoocP, the theoretical framework at 
finite temperature is less well developecP. Some representative studies in this regard include: a study of the effective 
potential for a weakly interacting Bose ga^, application of the background field method to study the interacting Bose 
gas at finite temperature and densitj^, a study of a two-species homogenous and dilute Bose gas using the effective 
potential^; and an investigation of Bose-Einstein condensation of Feshbach molecules at finite temperatures using 
a mean-field approaclP^ - a comprehensive review of various studies is also availabl^S] jj^ ^jjig paper we describe 
another approach to the study of a weakly interacting Bose gas at temperatures below the critical temperature. Based 
on a perturbative diagrammatic expansion of Green functions, the method invokes a self-consistent condition related 
to the minimization of the free energy. As will be seen, physically reasonable results are obtained in regions where a 
perturbative expansion is expected to be accurate. 

The paper is organized as follows. In Section |lT] we give a general overview of the procedure and estimate the 
thermodynamic parameters of the system under study. The algorithm itself and approximations made are discussed 
in Sections |IIB and II C followed by a discussion of the choice of the system parameters in Section II D Section III 



presents our findings and a discussion of the effects and features observed. In Section [rV] we give some conclusions. 

For reasons that will be discussed shortly, perturbative corrections in this paper are performed within the thermo- 
field dynamics (TFD) framework^^. Since this technique may not be as well known in this context, we discuss it in 
the Appendix in some detail: we first review its formal construction in Sections |A 1| and |A 2[ followed by a discussion 
of the Feynman diagrams in Section A 3 and a general statement of the problem of the weakly interacting Bose gas in 



Section [A 4[ The Appendix finishes with the formulation of the TFD Dyson-Beliaev equations, which serves as the 
starting point for considerations in the main body of the paper. 



II. METHOD 



A. Thermo-field dynamics formalism for a Bose-condensed gas 



In order to incorporate temperature effects, we will treat the system under study as an open one, exchanging 
energy and matter freely with a thermal reservoir, using the techniques of the field theory at finite temperature 
and densitjff^ti^. There are two main types of finite temperature formalisms: the imaginary-time, or Matsubara, 
formalism, and real-time formalisms. The path-integral formulation is a useful framework for seeing how these two 
formalisms are relatecP^. The imaginary-time formalism is straightforward, but an analytic continuation to real time 
must be done when the calculation of interest warrants it. The real-time formalism works directly in real time, but 
at the expense of a doubling of the number of degrees of freedom, which can make higher-order calculations tedious. 
The choice of which formalism to employ is, to a large extent, a matter of personal preference, as in the end both 
formalis ms mus t lead to the same physical results. We have found a real-time formalism, based on an operator-based 
approactPE^IiS, to be convenient for the lower-order calculations done here, as a straightforward interpretation of the 
various terms in the perturbative expansion is possible. 

The particular formalism we choose is the thermo-field dynamics (TFD) frameworlP^. Being a real-time formalism, 
a doubling of the degrees of freedom is needed at finite temperature, resulting in a 2 x 2 matrix propagator. The 
diagonal elements of this matrix are Green's functions which describe the propagation of excited atoms, while the 
off-diagonal elements are responsible for the exchange interactions. This propagator matrix is subject to the Dyson- 
Beliaev equations, expressed in terms of a 2 x 2 matrix self-energji22] These equations are general and coul d, in 
principle, be applied to any Bose-condensed system without reference to whether it is weakly interacting or nolp'^. 
However, except in special cases, a closed-form solution is generally not possible, and approximation schemes are 
often invoked. For a perturbative expansion, choosing an appropriate subset of Feynman diagrams contributing to 
the self-energies constitutes a starting point, which we describe in the next Section. 
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Procedure 



We start with the Dyson-Bchaev equations in TFD for the full propagator matrix, written in the following form 
(see Eqs. ( |A70||A72| ) in the Appendix): 



g(p) = g(°)(p) + g(°)(p).s(p).g(p), 



where 



G(o)"^(p) 







G(o)"/3(p) 



(1) 

(2) 
(3) 
(4) 



Here G and G'-^' are, respectively, the full and free TFD Greens functions (the latter is given by Eq. (A24|), while 
the elements of the E matrix are the self-energies. Due to the doubling of degrees of freedom, all of these functions 
have a 2 X 2 matrix structure, indicated by contravariant indices {a,/?}. Covariant indices {1,2} distinguish between 
the normal and anomalous types of functions, as described in the Appendix. The following identities help to reduce 
the number of independent components: 



(5) 



Eqs. ([!]) could be solved in principle with respect to the full propagators ([2|, expressed in term of the self-energy 
functions Q. Except in special cases, however, a closed form solution is typically not possible, in which case some 
type of approximation is often used. In this work, we take into account only the self-energy functions which correspond 
to the Bogoliubov approximation (see Fig.[T]). 



-<s> 



FIG. 1. First order corrections to the self energy functions of the interacting Bose gas in the Bogoliubov approximation. Solid 
lines correspond to propagators of the excited particles, wavy lines to the condensed phase ones. Dashed lines represents the 
interaction potential. Details of the notation is found in the Appendix. 



In this approximation one neglects all of the interactions between the excited atoms, with only those with the 
condensate taken into account. With this, we may solve ([l]) and find the propagators as functions of two unknown 
parameters: the condensed phase density no and the chemical potential fi. In TFD, one may express any physical 
parameter of interest via the G\\ matrix element of the propagator. In particular, the density n can be written as: 



n = no + 



J duj dp G"(cj,p;7io,/i) 



(6) 



which gives an implicit dependence of uq on fi. In order to completely characterize the system we also need, therefore, 
an independent way to determine the chemical potential itself. A common practice here is to use the Hugenholtz- 
Pines theorem^, which relates this parameter to the self-energy functions of the systerd^. However, we proceed in a 
different manner and require that the physical values of the condensate density and the chemical potential in a state 
of equilibrium should minimize the free energy F of the system. In the TFD formalism, interaction corrections to this 
potential are given by a set of vacuum diagrams in which all internal propagators are of G^ii'^^ type. As in the case 
of the self-energy, we maintain only the minimal possible subset of diagrams required, which is depicted in Fig. [2j 
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FIG. 2. First and second order corrections to the free energy of the interacting Bose gas in the Bogoliubov approximation. 



Thus, for example, in the first order of approximation the free energy reads: 

1 



where 



kT 
'A3 



55/2 (exp(//)) 



(7) 



(8) 



is the free energy of the ideal Bose gas. Here k is the Boltzmann constant, T is the temperature, U is the interaction 
potential, (75/2(2;) is the Bose-Einstein functiorP^l, and A — h/{2TTmkT)^/'^ is the thermal wavelength, which, in turn, 
depends on the mass m of an atom. 

In order to keep the procedure consistent, we take into account the interaction corrections to the self- and free 
energies which are of the same expansion order simultaneously. Our task thus consists of the following steps. 

• Given a subset of contributions to the self-energy, we iteratively solve ([6]) to find uq as a function of /i, T and n. 

• This value of tiq is then used to determine the free energy F. 

• The chemical potential is adjusted and this process repeated until the minimum of the free energy is found, 
which thus corresponds to the physical state of equilibrium. 



Once the chemical potential and the condensed phase density as functions of temperature are determined, the rest of 
the parameters of interest follow by applying standard thermodynamic relations. 



C. Discussion of the method 



At first glance, this approach might be seen to suffer from the drawback that the use the Bogoliubov self-energies 
to explore the thermal properties of the system in a strict perturbative expansion is only accurate in the vicinity of 
T = 0. However, the determination within our procedure of the chemical potential involves a partial resummation of 
terms in the following sense. Let Fn be the free energy in the nth-order of approximation, given by the sum of the 
zero order free energy Fq and corresponding loop corrections: 

F„(no, ^l) = Fo(no, ^l) + F^^\no, /i) + . . . + F^"\nn,fi). 

The chemical potential is determined by minimizing the entire right-hand side of this function, and so places terms of 
different loop order on the same footing. Thus, a partial mixing of loop order occurs, but the minimization of the free 
energy makes the outcome of this procedure thermodynamically consistent. The subset of the self-energy diagrams 
arising in the Bogoliubov approximation is thus regarded as a minimal starting point in this procedure, and we might 
anticipate this to potentially provide reasonable results at temperatures below the critical temperature. Of course, 
the ultimate test is to examine the results obtained via this procedure. 



D. Model parameters 



We choose liquid ^He as a representative system, and so take the atomic masses and molar volume to be m = 
6.65 X 10~^* g and V = 27.6 cm^/mole respectively^S. Since the momenta involved are sufficiently small due to the 
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low temperature of the system, we assume that the Fourier components of the pair-interaction potential L'^(p) could 
be replaced by the value J7o at p = (W. Thus, the form of the repulsive interaction potential we use is given by 

Uo = , 9 

m 

where a is the s-wave scattering length. In order to use Uq as a loop expansion parameter, the value of a is chosen 
such that 

?| - 0-02, (10) 

where Tc ~ 3.13 K is the critical temperature of the system. In other words, the interaction energy between the 
particles is assumed to be small compared to the thermal energy. Furthermore, unless otherwise stated, we will assume 
that the value of Uq is unity if it satisfies the assumption of (10 1. 



III. RESULTS AND DISCUSSIONS 



In this Section we present the results of calculations of the first and second order interaction corrections to ther- 
modynamic functions of the ideal Bose gas within this framework. 

In Fig. [3] we examine the dependency of the condensate density on temperature in the zeroth, first and second 
orders of approximation from ([6]). 

no/n 




12 3 



FIG. 3. The order parameter as a function of temperature. Red, blue and green lines correspond, respectively, to the zeroth, 
first and second orders of interaction corrections at unit interaction strength. 



Two features of interest are found here. The first is that the second order correction does not alter the overall 
picture too much compared to the first order one in the whole range of temperatures considered excluding, possibly, 
the region close to the critical. This is what one expects from higher-order loop corrections, and thus provides a 
degree of assurance that the perturbative expansion is valid in this region of temperatures. The second feature is that 
it is difficult to draw any definite conclusions about the system's properties in the critical region, since the procedure 
we have followed is not stable there. This is not a refiection of a particular calculational algorithm, but refiects the 
overall drawback of all theories that try to approach the critical region perturbatively by taking into account only a 
limited subset of self-energy diagrams into account-; all such approximations break down near the phase transition 
due to infrared divergencies. This problem might be cured by a selective resummation of diagrams from all orders 
of perturbation theory, but since we do not do that here, we cannot make any definite conclusions about the critical 
region, but rather limit our observations in this region to only general statements. 

Thus, we conclude that much of the information concerning the infiuence of the weak interaction on the order 
parameter in the whole range of temperatures excluding the critical region is contained already in the first order of 
approximation. Let us examine this in more detail. 
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In Fig. |4] the first order interaction correction noi from ([6| to the condensed phase density as a function of T/Tc 
is represented for three different values of the interaction strength. 



no i/n 




0.02 



FIG. 4. The first order interaction correction to the condensed phase density as a function of T/Tc- Green, red and blue lines 
correspond respectively to the interaction parameter Uo equal to 1/2, 1 and 2. 



First, we notice, that noi is positive everywhere except in a region infinitesimally close to T = 0. Thus, at very low 
temperatures the condensed phase density of the weakly interacting Bose gas is lower than that of the ideal gas, even 
at T = 0. This effect is a purely entropic nature and appears due to the additional enthalpic term contributed by the 
interaction. This term increases the entropy, which is physically equivalent to a decrease of the number of condensed 
particles. 

In a strict loop expansion, the region where npi stays negative determines the range of validity of the Bogoliubov 
approximation. However, as discussed in Section II C the partial mixing of loop orders makes possible potential 



conclusions about the system beyond the scope of the original Bogoliubov construction. Except in an infinitesimal 
region near T = 0, the condensed phase density of the interacting gas is higher than of its non-interacting counterpart. 
An explanation of this could be as follows. The physics of a non-ideal Bose gases is determined by the outcome of 
the competition between the entropy and enthalpy. While the first favours more particles in a normal phase by 
increasing the number of accessible microstates as the temperature increases, the second tries to put more atoms 
into the condensate by the use of repulsive interactions. The point where these two mechanisms balance each other 
corresponds to the state of thermodynamic equilibrium. The positivity of noi gives us a clue that the interaction term 
gives a greater contribution to the enthalpy than to the entropy in the whole range of temperatures, starting from 
the ones just outside of the Bogoliubov domain up to near T ^ Tq, where the density correction goes to zero. This 
seems to imply that the thermal effect is not as significant in an interacting Bose gas as in an ideal gas. 

As is clear from Fig. |4] the bigger the interaction strength, the bigger is the rioi correction. Qualitatively, this 
dependency is nonlinear, and two-fold. Firstly, as we observe in Fig.|4j an increase in the magnitude of the interaction 
is not accompanied by a proportional increase of the value of the correction itself. Secondly, the bigger the value of 
Uq, the greater its change impacts ngi. One may directly observe this effect from the relative distances between the 
curves while keeping T/Tq fixed. 

It is also worth noting that the maximum of Hqi shifts slightly to higher values of temperatures as L^o increases. 
This suggests that the enthalpy is able to dominate the entropy longer the stronger the interaction is between the 
atoms. 

What is also interesting is the absence of a symmetry between the curves with respect to the midpoint T/Tc = 1/2: 
the bigger the interaction strength, the bigger the asymmetry. As discussed above, this is due to the nonlinear 
competition between the entropy and enthalpy, which leads to different slopes on the left and right sides of the 
maximum. 

The plot of the dependency of tiqi on T/Tc functionally and on Uq parametrically captures most of the physics. 
We can also extract some limited information about the behaviour of the system as the temperature approaches the 
critical temperature by examining the first order interaction correction to the specific heat Cyi per particle as a 
function of T/Tc- This is shown in Fig. [sj 
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Cy 1 

N k 




FIG. 5. The first order interaction correction to the specific heat per particle as a function of T /Tc- Green, red and blue lines 
correspond respectively to the interaction parameter f/o equal to 1/2, 1 and 2. 

One readily observes the sharp rise of the specific heat as the temperature nears the critical temperature, with 
the bigger the interaction parameter, the sharper the rise. This is the expected qualitative behaviour of interacting 
systems near the phase transition point. 

One final feature of this approach should be noted. It turns out that in the first order of approximation the chemical 
potential vanishes. This is in contrast to resuhs obtained within some other approaches where the chemical potential 
of the ideal Bose gas acquires a small positive correction due to the interaction at first order. In our approach, the 
second-order correction appears significant in this case, as with this, the chemical potential acquires a small negative 
value, comparable in order of magnitude with the interaction strength. The magnitude is maximal at T = and 
quickly drops to zero in absolute value as the temperature rises, as seen in Fig. [6j 

k Tc 
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-0.02 
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-0.04 



0.2 0.4 0.6 0.8 1 



FIG. 6. The chemical potential in the second order approximation with unit interaction strength. 

This feature of second-order effects being significant when first-order effects vanish is not rare; this is found, for 
example, in the hard thermal loop resummation of Braatcn and PisarksP^and in certain effective potential calculations 
involving electroweak interaction^2ll_ 

IV. CONCLUSIONS 

Based on a combination of a loop expansion of the self-energy and a self-consistent minimization of the free 
energy, we have shown how thermal properties of the weakly interacting Bose gas could be explored within a field 
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theoretic framework. Despite the fact that our procedure uses only the minimal subset of self-energy and free energy 
contributions, this allows one to evaluate the major effects caused by the interaction already in the leading loop 
expansion order. A notable feature found is an increase in the condensate density, compared to the non-interacting 
case, in the whole temperature range excluding a small region near T — 0, and a sharp rise of the specific heat near 
the phase transition point. These effects clearly demonstrate that the thermal properties of ideal Bose gases could be 
changed significantly by the interaction. 

However, as with all such perturbative expansions, the physically interesting region near the critical temperature 
remains inaccessible by our method. Summation to all orders of a broad class of self-energy diagrams is likely required 
to make reliable conclusions about phenomena taking place near the phase transition pointP. 



Appendix A: The TFD framework 



We have chosen the TFD approach as a framework to perform the calculations in this paper due to the physical 
interpretation possible in a perturbative diagrammatic expansion. Although a doubling of the degrees of freedom is 
entailed, for lower orders this complication is not extreme, and in any event, physical quantities should be independent 
of the formalism used to calculate them. 



1. Formal construction of TFD 

The TFD formalism is based on the following argument^^^l. Let us imagine an open system, freely exchanging 
particles and energy with a heat reservoir, maintained at constant inverse temperature /3. The process of annihilation 
of a quanta of energy w(k, /3) in such a system can be regarded as the action of an operator a(k) which satisfies the 
Bose commutation relation: 

[a(k),at(i)] ^j(k-l). (Al) 

However, unlike the zero temperature case, the operator a(k) has an internal structure which incorporates two 
independent kinds of thermal effects: 

• annihilation of a quanta in the system; 

• creation of the (Dirac) "holes" maintained by the reservoir; 

When the latter process takes place, we say that a a^-quantum (hole) with negative energy — ?iaj(k) and negative 
momentum — ?ik is created. The former process will be described by the annihilation operator a(k, Since a(k) 
now performs two independent operations, it can be represented by the linear combination: 

a(k) - c(k, /3)a{k, /3) + d{k, P)a^ (k, (A2) 

where c(k,/3) and d(k, /3) are certain c-numbcr functions. Similarly, the creation operator is 

a^(k) = c*(k,/3)a^ (k, /?)-!- d*(k,/3)5(k,/3), (A3) 

By construction, the newly introduced temperature dependent operators satisfy Bose commutation relations. This 
suggests that relations ( A2|A3 ) could be regarded as a canonical transformation, which implies 

c(k,/3)c*(k,/3) - d(k,/3)d*(k,/?) = 1. (A4) 

The phase factors of the above coefficients could be absorbed into those of the respective operators, so we can choose 
c(k, /3) and (i(k, /3) to be real. 

In order to establish their explicit form, let us use the following considerations. The operator Q;^^(k) takes care of 
the excitations in the system. Therefore, the quantum number density operator is 

Af(k) = a^(k)a(k). (A5) 

Since the model we consider is in the grand canonical ensemble, we may require that the vacuum average of this 
operator should be equal to the corresponding value of the number density in statistical mechanics: 

n(k) = iM!(o,/3|^(k)|0,/3) = /s(w), (A6) 
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with 



1 

g/3u;(k,/3) _ l' 



This gives 



(A7) 



(A8) 



Eq. ( |A6[ ) constitutes the axiom which determines the temperature of the system in a thermal equihbrium state. It 
can easily be extended to any operator consisting of a and a\ say A{a,a^). Thus, the vacuum expectation value 
{0, f3\A{a,a^)\0, P) is equal to the ensemble average of A{a,a^). In this way, the statistical mechanics of a quantum 
many-body system becomes a quantum field theory realized in a temperature-dependent Fock space. This formalism 
is called thermo-field dynamics. 

The form of Eqs. (A2 A3) suggests that tilde and non-tilde operators are related to each other. Indeed, there 



is a general TFD "tilde-conjugation" operation which is valid not only for the elementary creation and annihilation 
operators, but also for arbitrary compound functions Oi and O2: 



O1O2 = 0i02, 

CiOi~C2C'2 =c{di+ c;02, 

6^0. 

Constructing the free TFD Hamiltonian Hq is now straightforward: 

Ho ^ Ho — Hq. 

with 

Ho = hJ d k w(k)a^(k,/3)a(k,/3). Ho = h J d k w(k)5^(k, /3)5(k, /3) 
or, equivalently, with the use of the canonical transformation 

Ho = n J dk w(k)a^(k)a(k), Ho^fi j dk w(k)(5^(k)<5(k). 



(A9) 
(AlO) 
(All) 

(A12) 
(A13) 

(AM) 



2. Green's functions in TFD 



Let us now consider the simple example of a free field at finite temperature. Suppose such a field (j)^ satisfies the 
equation \{d)4f' = 0, with 



A(a) = z|-^v). 

The expressions for the tilde- and non-tilde Haniiltonians read: 

Ho^ J dx 0°^w(V)0°, 
i?o = [ dx 0°^a;(-V)0°, 



(A15) 

(A16) 
(A17) 



where the relation i^*(V) = w(— V) from the tilde-conjugation rules ( A9|All ) has been used. The corresponding 
equations of motion 



d 



i;i-0°(x) = [0"(x),iJo], 



(AI8) 



6"(x),i/o 



(A19) 
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have the formal solution: 



0°(x) 



Vn 



(27r)3/2 



J dk Q;(k) exp[ik • x — ioj{\i)t], 



dk a(k) cxp[— ik • x + iaj(k)t]. 



(27r)3/2 

It is useful to simplify the notation by introducing the following thermal doublet symbol: 



(A20) 
(A21) 

(A22) 



The causal two-point function is then defined by 



= e{t, - ty){o, i3\r (x^^^ |o, /3) + e{ty - t,){o, /3|0^t (y)</,"(x) |o, /?) 



(A23) 



where the symbol T denotes ordering along the TFD time contour in the complex plan^^, and the average is over 
the temperature-dependent vacuum. Using (A20 A21| and the canonical transformations (A2 A3), we find that in 



the momentum representation the above function reads: 



G(°)°'5(fco,k) = 



/ C%{uj) CB iuj)dB (^) _ CB{uj)dB{L0) \ 

fco — w + ie kg — uj — ie ko — uj + ie kg — lu — ie 



CB{i^)dB{uj) _ cb iuj)dB juj) 
\ ko — U} + ie ko — u! — ie ko — UJ + ie ko — uj — ie 



dlic.) 



c%iuj) 



(A24) 



where lu — a;(k) and cb(w), ^^(a;) are given by (A6). This free propagator (A24) is the building block for the 
development of perturbation theory, which we now discuss. 



3. Feynman diagrams in TFD 

Perturbative calculations in TFD are similar to those at zero temperature, and indeed, one of the virtues of the 
formalism is that the machinery of renormalization and the renormalization group equations are readily incorporated. 
The doubling of degrees of freedom discussed above leads to a 2 x 2 matrix structure of propagators and self energies. 
This feature implies that in higher orders a much larger number of diagrams has to be taken into account, compared 
to the vacuum theory, and calculations may become cumbersome. However, for calculations to low orders, this is not 
a serious problem. 

In general, Feynman diagrams turn out to be topologically and combinatorically identical to those of the corre- 
sponding zero-temperature quantum field theory, but generalized to a two-component formalism. The first (type 1) 
component represents, by construction, the "physical" field, describing propagation of excitations in the system. The 
second (type 2) component corresponds to a "thermal ghost field" , corresponding to holes, maintained by the reser- 
voir. In TFD, two types of vertices occur: one type describing the interactions of the type 1 fields, and the other of 
the type 2 fields. The Feynman rules for these two types of fields differ only by a sign in the case of even interactions. 
There is no direct coupling between the two types of fields, but they can propagate into each other because of the 
non-diagonal 1-2 and 2-1 elements of the propagator matrix. The 1-1 and 2-2 components correspond to physical 
and ghost propagators. Therefore, external lines of physical Green's function are always of type 1. In order to find 
a particular n-point Green's function in momentum space, one must draw all topologically distinct diagrams with n 
external points of type 1, and then sum over internal vertices of types 1 and 2. The ultimate goal is to compute the 
1-1 matrix component of the Green's function, which can then be used to calculate physical quantities of interest. 



4. Weakly interacting Bose gas in TFD 

Let us now formulate the general problem of a weakly interacting Bose gas within the TFD framework. 

The system we study is a dilute gas consisting of N atoms obeying Bose statistics, enclosed in a box of volume V and 
interacting through a two-body potential C/(x — y). The interaction is assumed to be weak which, together with the 
low density of the gas, allows us to neglect higher-order interactions and treat the system's properties perturbatively. 
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For simplicity, the atoms are assumed to have zero net spin, and are described by the boson field operators ?/'(x), ■0^(x) 
subject to the Bose commutation relations: 



[V'(x),V(y)] -0, 

[V'^(x),V;^(y)] -0, 

[^(X),^t(y)]^5(^_y) 



(A25) 
(A26) 
(A27) 



The system is immersed in a large heat reservoir, with which it can freely exchange both energy and particles. In a 
state of equilibrium both the system and the reservoir possess the same chemical potential /x and inverse temperature 
/?. In TFD it is assumed that the reservoir, unlike the system, maintains the negative-energy holes, represented by 
the operators obtained by the tilde-conjugation operation of ( A9||A11 ). They satisfy the same commutation relations 
as the regular operators, while all commutators between the tilde- and non-tilde operators vanishing by construction. 
In equilibrium, the system contains a mixture of particles and holes. Therefore, its Hamiltonian H is given by 



H = K -K, 



where 



K = Hq- fiN 



2m 



(A28) 



(A29) 



is the regular zero-temperature Hamiltonian describing the particles of the system and K is its tilde-conjugate, 
describing the holes: 



K ^ Ho- nN ^ J dx ^^x) 



^(x) 



2rn 

dyV;t(x)V;^(y)C/(x-y)V3(y)V3(x). 



(A30) 



The chemical potential /i is chosen such that the vacuum average of the atomic number operator is equal to the total 
number of particles in the system: 



(N) = N. 



(A31) 



With the Hamiltonian H, one may introduce the time-dependent Hcisenberg picture for any Schrodinger operator 
0(x) and its tilde-conjugate 0(x): 



0(x,i) = e*^"'*0(x)e-*^'*, 
0(x,i) = e-*-^*0(x)e*-^*. 

In particular, the field operators have a time-dependence given by 

^(x,i) = e^^V(x)e^^^*, 

iKt. 

kt 

iKt 



V'^(x,i) = e-'^*tlj\x)e 
ij{x,t) = e-''^'^{x)e 
ij\x,t) ^ e:^*^p\x)e- 



In order to simplify the notation, we introduce the thermal doublet symbol (see (A22)) 

^i(x,i) \ _ f V(x,<) 



The single-particle Green's function is then defined as: 

G"^(x,i;x',t') = {0,m^"i^,t),i^^H^',t')]\0,/3), 



(A32) 
(ASS) 

(AS4) 
(ASS) 
(AS6) 
(AS7) 



(ASS) 



(AS9) 
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For our purposes, it is more convenient to work in the single-particle momentum representation. In this represen- 
tation, the field operators ip{yi),tp{x) can be written as 

V'(x) = . / dk ak(Ak(x), (A40) 



^(x) = J dk ak<^k(x), (A41) 

where ak(a|j.) and ak(ak) are the creation (annihilation) operators for a particle and hole in single-particle states 
(/)k(x) and ^k(x), respectively. These obey Bose commutation relations: 

[ttk, flk'] = 0, [aw, al,] = 5k,k'; (A42) 

[ttk, flk'] = 0, [ok, dl,] = (5k,k', (A43) 



with all commutators between the tilde and non-tilde operators vanishing. 

Just as in the case of fields, we may introduce the thermal doublet notation for the operators Ok and Ok and then 
write the Green's function G"^(x, f;x',t') in the single-particle momentum representation as 

G"'5(x,t;x',t') = ^^JdkJ dk'4>Zi:^),^iH^'mP\TK{t),ai}{t')]\0,P}. (A44) 

In the case of a uniform system governed by a time-independent Hamiltonian, the single-particle states are plane 
waves 

(^k(x) = exp(ik • r), (A45) 



V 

0k (x) = exp(-zk • r). (A46) 



Similarly, G"^{x, t; x', t') now depends only on the differences x — x' and t — t', so that 

G"'^(x,t) =y"dk e*-^G"'^(k,<), (A47) 

where 

G«^(k,i) = (0,/3|TK(t),a^t(o)]o,^) (A48) 

is the single-particle Green's function in the momentum representation. 

The remaining task is to calculate the Green's function G"^{k,t). For a system of Bose particles, this task is 
complicated by the possible phase transition to a Bose-condensed state by spontaneous symmetry breaking below a 
certain critical temperature Tq- This is signalled by a non- vanishing average 

(0,/3|V"(x)|0,/3) = <i>^(x)^0, (A49) 

The function $o (x) is often referred to as the macroscopic wave function of the condensate and in general has an 
amplitude and phase. In a uniform system and in the absence of any supercurrent, we can take $o (x) to be real and 
independent of position. In this case, $0 (x) is equal to the square root of the condensate density tiq. We are then 
led to separate the boson field operator ip"{x.) into two parts: 

V'" (x) EE a^/y + V!" (x) , (A50) 
^/'"^(x) = a'^yv^^^ + V>^(x). (A51) 

The commutator of Uq and a^^ is unity, while their product is of order A''. Therefore, we may neglect the operator 
nature of o" and a^^ and replace them by the c-number ^/Nq. This procedure, known as the Bogoliubov prescription, 
is appropriate when the number of particles in the zero-momentum state is a finite fraction of N. The error thus 
introduced is of the order 0{V~^) and vanishes in the thermodynamic limit. The new field operators ^"(x), ^"'''(x) 
describe the normal phase atoms, satisfy Bose commutation relations, and have vanishing vacuum averages: 

(0,/3|Vi"(x)|0,/3) = (0,/?|^"t(x)|o,/3) = 0. (A52) 
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With this, the thermal Green's function now becomes: 

G"'^(x,t;x',0 - -no + G"'^(x,t;x',t'), 
where G"'^(x, x', i') is its non-condensate part defined by 

G"^(x,t;x',t') ^ -{Q,mrM,i^P\^,t')]%P). 
In the momentum representation for k 7^ 0, we have 

G"^(k,i)--(0,/3|rK(t),a^^(0)]|0,/3). 



(A53) 



(A54) 



(A55) 



The Bogohubov prescription modifies the Hamihonian (A28) in a fundamental way. In order to simplify matters, we 



will concentrate on the zero-temperature part K for the moment, making a straightforward generalization to TFD 
later on. After all of the above transformations, K reads 



I 7 



(A56) 



where the interaction term V can be separated into the following eight distinct parts: 



k 

V2 = ^fio ^ L/(k)aj^a_k, 

k 

V3 = ?^oX!^('*)"k"k, 

k 



^^6 = 



k 

1/2 

072 X! t^(q)ak+qakaq, 

k,q 

1/2 

172 51 t^(q)akaqak+q, 

k,q 



^7 = ^ H C/(q)ak+qaL_qak'ak 



k,k',q 



with 



t/(q) = y e'i-'^(7(x). 



(A57) 
(A58) 

(A59) 
(A60) 
(A61) 

(A62) 

(A63) 
(A64) 

(A65) 



In a normal system (no — 0), only V-j is present. We also note that K has no term containing a single or Ok 
because these would violate momentum conservation. This is in agreement with (A52), which gives in a momentum 
representation 



(0,/3|ak|0,/3) = (0,/3|aj,|0,/3)-0, (k ^ 0). 



(A66) 



The fact that there exist interaction terms such as Vi implies the non-conservation of particles in the system 
described by non-condensate operators ak and due to exchanges between non-condensate atoms and the condensate 
atoms. This changes the way the proper self-energies are defined. There are three distinct types of them now, indicated 
in Fig. [7} 
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1,1 (P) Ii2(P) ^21 (P) 



FIG. 7. Proper self-energies for a Bose-condensed system. Straight lines correspond to the normal phase particles, with wavy 
lines associated with the condensed phase. 



For Sii there is one particle line coming in and one coming out. The others components have two particle Unes either 
coming out (S12) or going in (S21), and reflect the new features associated with the existence of a Bose condensate. 
We then introduce two new Green's functions: 

Gi2(k,<) = -(0,/3|Ta_k(t)ak(0)|0,/3), (A67) 

G2i(k,t) = -(0,/3|Ta|,(t)aLk(0)|0,/3). (A68) 

The functions G12 and G21 are usually called "anomalous" Green's functions, representing the disappearance and ap- 
pearance of two non-condensate particles, respectively. The normal Green's function is denoted as Gn and represents 
the propagation of a single particle 

Gii(k,t) = -{0,l3\Takit)al{0)\0,l3). (A69) 

These three Green's functions are shown in Fig.|8] where the arrows indicate the direction of momentum of the atoms 
involved. 

P -p p p -p 

Gii(P) Gi2(p) G2i(p) 



FIG. 8. Green's functions for a Bose-condensed system. 

The Dyson equations connecting the above free energies and propagators is found in standard textbooks. For the 
case of zero-temperatur^^, they where first written down by Beliaev, and are illustrated diagrammatically in Fig. [9] 




FIG. 9. The Dyson-Beliaev equations for a Bose-condensed system. Thin and thick straight lines denote, respectively, the free 
and full propagators of the normal phase particles. 
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The generalization of this to TFD is straightforward and in the momentum representation reads: 

G«f (p) = G(°)"^(p) + G?7(p)S7f (p)G(")^^(p) + G^7(p)S72'(p)G(°)^^(p), (A70) 

= G77(-p)S72*(p)G(°)^^(p) + G72^(p)E7f (p)G(°)^''(p), (A71) 

G^f (P) = G77(p)Elf(p)G(°)^^(-p) + G^7(p)E^2'(~P)G(°)'^(p). (A72) 

Here the TFD indices a, /3 take the values of 1 and 2 and summation over repeated indices is implied. The n otation 



P = {Po,p) is employed and the expression for the unperturbed Green's function G'^^^°'^{p) is given by (A24|. There 
are some useful relations between the anomalous proper self-energies and propagators: 

G72''(p) = G^f(-p), G?f(|,) = G7f(-p), S^^^b) = S^f I]^2^b) = E7f(-p) (A73) 



Eqs. ( A70|A72 1 could be solved with respect to the propagators in terms of the exact proper self-energies and unper- 



turbed Green's functions. Therefore, they are completely general and can be applied to any uniform Bose-condensed 
fluid, liquid or gas, without any reference to whether the underlying system is weakly interacting or not. How- 
ever, except in special cases, a closed-form solution is generally not possible, which leads to consideration of various 
approximation schemes. Many such schemes, including ours, use Eqs. ( A70|A72 ) as a starting point. 
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